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probability distribution law one can determine the threshold U x to which the value is to be 
compared to result in a relatively small acceptable probability of false alarm. For example, an 
appropriate threshold is approximately P fa =0.01 . . .0.05, where X 2 is calculated by equation (11). 
If the threshold is exceeded, then it is determined that an anomalous measurement has been 
detected. Once the anomalous measurement is detected, a determination of the particular 
channel having the anomalous measurement is made by repeatedly removing channels from the 
X 2 calculation in a manner similar that as described above in accordance with the first described 
embodiment of the residual analyzer 320. 

A third embodiment of the residual analyzer is also based on criteria ^ , but locates the 
channel with the anomalous measurement in a different way. Once the threshold has been 
exceeded, it is assumed that an anomalous measurement exists in one channel only and the 
channels are eliminated one by one. However, upon eliminating a particular channel and 
recalculating X* > if the threshold remains exceeded using the remaining n=N-l channels, then the 
particular channel is returned to the set of channels used in the calculation of x 2 and a second 
channel is eliminated. Again, is recalculated with the remaining n=N-l channels. A first 
search continues in this manner until either the one channel having the anomalous measurement 
is found, or after eliminating all channels and the threshold continuing to be exceeded, it is 
determined that more than one channel has an anomalous measurement. In the latter case, a 
second search is initiated now eliminating channel pairs and computing X* using n=N-2 
channels. A second search continues in this manner until either the two channels having the 
anomalous measurements are found, or after eliminating all channel pairs and the threshold 
continuing to be exceeded, it is determined that more than two channels have anomalous 
measurements. In the latter case, a third search is initiated now eliminating channels in groups 
of three and computing x 2 using n=N-3, and searching again as described above. Searching 
continues in this manner until either all the channels that contain anomalous measurements are 
found, or the number of the channels remaining is equal to two. 

The three above described embodiments of the residuals analyzer 320 are appropriate 
where the integrated discriminator 308 generates the average estimate of phase mismatches A* 
according to equation (3). The operation of the residual analyzer 320 is altered as follows if the 
integrated discriminator 308 generates the average estimate of phase mismatches A* according 
to equation (5). First, the search is stopped when the number of remaining channels is five. 
Second, the degree of freedom for the calculation of the threshold in the second and third 
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residuals analyzer embodiments described above is n-4. Third, each time a new channel or 
channels is eliminated from the calculation, the integrated discriminator 308 recalculates matrix 
G using the remaining channels. In yet another alternative, in the case in which the. receiver is a 
dual-frequency receiver and phase difference is calculated using equation (1), if the search stops 
5 when the number of channels reaches five, it may be continued by calculating phase difference 
with equation (2) rather than equation (1). This results in locating more channels with anomalous 
measurements. 

Following the analysis of the residuals, the decision module 322 makes a decision based 
on the analysis. In one embodiment, the decision may be to merely transfer to the navigation 

10 location determining module 206 (Fig. 2) the information so far obtained about the anomalous 
measurements so that the navigation location determining module 206 may consider this 
information in solving the navigation task. If the combined phase difference is defined by 
equation (1), then the anomaly indicator module can also calculate an estimate of cycle slips to 
allow the navigation location determining module 206 to correct for such cycle slip. The 

15 estimate * c s (i) is based upon rounding the residuals and may be calculated as follows: 
*c s (i) =0.5 round(5(i)/ 0.5), (12) 
where * c s (i) is the vector of estimates (in semi-cycles) of anomalous errors (cycle slips), 
5(i) is the vector of residuals for the deleted channels, and round (...) is the operation of 
rounding to the nearest integer. 

20 In addition, the anomaly module can transfer to the navigation location determining 

module 206 information about channel indicator alarms 324 that mark the appearance of strong 
perturbations in the channels. These alarms are received from the channel indicators which 
operate as described below. 

One skilled in the art would recognize that the navigation location determining module 

25 206 may use the information received from the anomaly indicator module 210 in various ways. 
For example, the navigation location determining module 206 may miss ignore an i~th 
measurement, eliminate corrupted channels or give them less weight, apply corrections to slips,* 
or use the measurements of all the channels in coordinate determination. 

We now describe the operation of the channel indicators. The task of the channel 

30 indicators is to detect disturbance in the satellite channels that may cause anomalous phase 
measurements and to mark (i.e., flag) the channels with an alarm. 



Page 19 of 39 



Docket No. 101 0-0001 -USA 



S(i) = A 6(i-l) + AS(i), (17) 
where A is a coefficient that determines time constant and the gain of the filter-integrator, where 
A < 1 . The recommended value for A is in the range A=0.995. . .0.999. In normal mode the 
integrated residuals represent small noise fluctuations. When an anomalous measurement 
5 appears the integrated residual of the corresponding channel increases and the sum of squares 
residuals between channels also increases. 

The residuals analyzer 516 and the decision module 518 may operate in a manner similar 
to that described above in accordance with the residuals analyzer 320 and decision module 322 
of Fig. 3. However, in the embodiment under present discussion, the operations of the residuals 
10 analyzer 516 and decision module 518 must be carried out with the integrated residuals. 

It is also possible to use the x* criteria (as described above) only to determine the 
existence of anomalies, and to use channel indicator alarms 520 to determine the location of the 
channel on which the anomaly has occurred. (It is noted that this alternative is also applicable to 
the first and third embodiments of the anomaly indicator.) In accordance with this alternative, all 
1 5 channels marked with alarm flags are eliminated from the anomaly indicator module's 
calculations as described above, and then the above described criteria is utilized for the 
remaining channels. Depending on the result, the following decisions can be made. 

If the number of channels remaining is not less than four and x 2 is less than the threshold 
Tl x , then in addition to transferring information about the corrupted channels to the navigation 

20 location determining module 206, the decision module 518 can also measure anomalous errors in 
those channels using equation (12) to apply cycle-slip-corrections to * c(i). After applying the 
corrections, the residuals analyzer 516 can again use the x 2 criteria to assess the corrections 
applied. Further, it is possible to continue verifying the validity of the channel alarms by, for 
example, returning certain channels previously eliminated from the anomaly indicator module's 

25 calculations, when the flags/alarms of the amplitude and angle indicators were different. It is 
noted that if it is impossible to satisfy x 2 criteria then i-th measurement should be considered as 
uncertain. 

Fig. 6 shows a functional block diagram of a third embodiment of the anomaly indicator 
in which initial ambiguities are deleted due to processing the increments of combined phase 
30 differences over the interval from initial to current measurements. In addition, this embodiment 
of the indicator module is based upon the readings of the channel indicators. As represented in 
block 602 the inputs of the anomaly indicator module are single differences of full phases 
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measured at the Base and Rover for each satellite. The combined phase difference generator 604 
generates the vector $ c(i) of combined phase difference as described above in accordance with 
equation 13. Operator 606 generates the increment vector as follows: 

A*c(i)=*c(i)- *c(0), (18) 
where *c(0) is the value of the combined phase difference at the initial moment (i=0), the initial 
moment being reset periodically by the initial time moments reset unit 608, with a period equal 
to approximately 50 . . . 100 seconds and also resetting upon the appearance of a new satellite 
signal. 

Slip corrections * c s (i-l) generated by the correction unit 616 (after a delay of one unit 
imposed by delay unit 618) are applied to the increments by operator 610 to generate corrected 
increments A* c (i) as follows: 

A$c (i)= A$c(i) - *c s (i-l) (19) 
Next, the channel indicator alarms 624 are analyzed by the channel indicator analyzer 622 and 
the channels marked with alarms are deleted from further processing by the integrated converter 
612. If the number of the channels remaining for processing is at least four, then the corrected 
increments A* c (i) of the remaining channels are provided to the integrated converter 612 to 

A 

generate A^ (i) as follows: 

A 

A<I> (i) = H(i) ■ G(i) * A* c c (i) (20) 
Here, matrices H and G (as described above) are calculated only for the remaining channels. 
Next, the vector of residuals 6(i) for the corrected increments is generated by operator 614 as 
follows: 

S(i) = A*c c (i)- A<*> (i) (21) 

The residuals are then provided to the correction unit 616 where cycle slip-corrections are 
generated as described above in accordance with equation (12). As described above, these 
corrections, after the imposition of a delay by delay unit 618, will be used for the calculation 
associated with the next measurement period. As seen, the correction loop in this embodiment of 
the anomaly indicator module is closed. In fact, the correction loop will work if alarms of the 
channel indicators appear. 

The decision module 620 operates in manner similar to that described above in 
accordance with the decision module 322 of Fig. 3. For example, the results obtained in the 
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Claims: 



1 1 . A method for detecting anomalous phase measurements in a satellite differential 

2 navigation system in which a first satellite receiver and a second satellite receiver each compute 

3 phase measurements on a plurality of satellite channels, said method comprising the steps, 

4 performed for each of a plurality of iterations, of: 

5 a) generating a combined phase difference vector; 

6 b) generating a phase mismatch vector representing the difference between said 

7 combined phase difference vector and an estimated combined phase difference vector; 

8 c) generating an averaged estimate vector representing the averaged estimate of said 

9 phase mismatch vector over said channels; 

10 d) generating a residual vector representing the difference between said phase mismatch 

1 1 vector and said averaged estimate vector; 

12 e) generating a vector of controlling signals by linear transformation of said residual 

1 3 vector and said averaged estimate vector; 

14 f) generating said estimated combined phase difference vector by successively storing 

1 5 components of said vector of controlling signals for each of said channels; and 

16 g) detecting anomalous phase measurements by analyzing said residual vector. 

1 2. The method of claim 1 wherein said combined phase difference vector * c comprises 

2 phase differences (<£c/..<i>c / ..<£cyv) for a j-th satellite channel, and wherein step a) of generating a 

3 combined phase difference vector further comprises the step of: 

4 calculating said phase difference vector * c as *c = * B -* R where vector * B 

5 comprises the full phases measured by one of said satellite receivers for eachj-th satellite 

6 channel ($ bj- . * bj- $bn) and vector * R comprises the full phases measured by the other satellite 

7 receiver for each j-th satellite channel rj. rn)- 

1 3. The method of claim 1 wherein said first and second satellite receivers are dual 

2 frequency (fl and f2) receivers and wherein said combined phase difference vector $ c 

3 comprises phase differences ($ci..$Cj..$c N ) for a j-th satellite channel, and wherein step a) of 

4 generating a combined phase difference vector further comprises the step of: 

5 calculating said phase differences $9 of said phase difference vector * c as 
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6 *c r (<t» B 2j - <t> R 2j ) - (*% - $%)(£> / fl) 

7 where 3> B 2 j and $ B ij represent the full phase measured by one of the satellite receivers at each j- 

8 th satellite channel at frequencies f2 and fl respectively, and 

9 where $ R 2 j and 3>% represent the full phase measured by the other of the satellite receivers at 
10 each j-th satellite channel at frequencies fl and fl respectively. 



1 4. The method of claim 1 wherein said step c) of generating an averaged estimate vector 

2 &ty further comprises the step of: 

3 calculating said average estimate vector A^ as 

4 N N 

5 A* = ( Ewj- A*cj)/( Ewj) 

6 j=1 j=1 

7 where A$Cj represents the components of said phase mismatch vector for a j-th satellite channel, 

8 Wj represents a weight coefficient for each j-th satellite channel, and N represents the number of 

9 satellite channels. 

1 5. The method of claim 1 wherein step c) of generating an averaged estimate vector A* 

2 further comprises the step of: 

3 calculating said average estimate vector A^ as A^= H ' G * A* c where H is a matrix of 

4 directional cosines, AS> c represents said phase mismatch vector, and G is a matrix defined by 

5 G=(H T R" 1 H)" 1 H T R" 1 where R is a covariance matrix of phase mismatches. 

1 6. The method of claim 1 wherein said step e) of generating a vector of controlling 

2 signals further comprises the step of: 

3 generating said vector of controlling signals in combined loop filters of each of said 

4 plurality of channels. 

1 7. The method of claim 6 where each of said combined loop filters generates a control 

2 signal ( U c j) for a j.-th channel according to the following: 

3 i 

4 UuCi^lTSjOc) 
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5 



1 



1 



th satellite channel BJ ..*bn) and vector * R comprises the full phases measured by the 



6 other satellite receiver for each j-th satellite channel ri..$ Rj ..$> m). 



23. The apparatus of claim 21 wherein said first and second satellite receivers are dual 
frequency (fl and f2) receivers and wherein said combined phase difference vector 4> c 
comprises phase differences (*c,..*Cj..*c N ) for a j-th satellite channel, and wherein said phase 

4 difference generator generates said phase differences * Cj of said phase difference vector * c as 

5 * Cj = (* B 2j - 4» R 2j ) - (*V* R ijXf2/fl) 

6 where 4> B 2j and * B ,j represent the full phase measured by one of the satellite receivers at eachy- 

7 th satellite channel at frequencies f2 and fl respectively, and 

8 where * R 2j and represent the full phase measured by the other of the satellite receivers at 

9 each>fA satellite channel at frequencies f2 and fl respectively. 

1 24. The apparatus of claim 21 wherein said integrated discriminator generates said 

2 averaged estimate vector as 

3 n N 

4 A* = ( Ewj-A*Cj)/( Iwj) 

5 j =l j=' 

6 where A* Cj represents the components of said phase mismatch vector for a j-th satellite channel, 

7 Wj represents a weight coefficient for eachy-//» satellite channel, and N represents the number of 

8 satellite channels. 



1 .25. The apparatus of claim 21 wherein said integrated discriminator generates said 

2 averaged estimate, vector A* as A*= H G A* c where H is a matrix of directional cosines, 

3 A* c represents said phase mismatch vector, and G is a matrix defined by 

4 G=(H T R" , H)" 1 H T R" 1 where R is a covariance matrix of phase mismatches. 



26. The apparatus of claim 21 wherein said joint loop filter generates said vector of 



2 controlling signals in combined loop filters of each of said plurality of channels. 



27. The apparatus of claim 26 wherein each of said at least one combined loop filters 



2 generates a control signal ( U cj ) for a j-th channel according to the following: 
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